Kink-antikink collisions in the cf)^ model 
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We study kink-antikink collisions in the one-dimensional non-integrable scalar (j)^ model. Although 
the single-kink solutions for this model do not possess an internal vibrational mode, our simulations 
reveal a resonant scattering structure, thereby providing a counterexample to the standard belief 
that existence of such a mode is a necessary condition for multi-bounce resonances in general kink- 
antikink collisions. We investigate the two-bounce windows in detail, and present evidence that 
this structure is caused by existence of bound states in the spectrum of small oscillations about a 
combined kink-antikink configuration. 

PACS numbers: ll.lO.Lm, 11.27.+d 



Introduction. Solitary wave solutions in non-linear field 
theories have been studied for many decades. Perhaps 
the simplest examples are the kink solutions which ap- 
pear in models in one spatial dimension with a poten- 
tial with two or more degenerate minima. These models 
have applications in condensed matter physics [H, field 
theory [3, Q and as simplified models of domain wall 
collisions in cosmology (^-[gI . The double well case is ex- 
emplified by the nonintegrable 0^ model. A particularly 
striking feature of this theory is the intricate structure 
of kink-antikink scattering, first observed numerically in 
the 1970s and much studied thereafter [sj-fist. For 
initial velocities above a critical value Vc ~ 0.2598 the 
two incident waves always escape to infinity after colli- 
sion, with the emission of some radiation. Below Vc, the 
incident waves generically become trapped, but there is 
also a complicated pattern of narrow 'resonance windows' 
within which they are again able to escape to infinity. 

The accepted explanation of the appearance of these 
windows, both in this model and in similar theories such 
as the parametrically modified sine-Gordon model [lo| . 
is that they are related to a reversible exchange of energy 
between the translational and vibrational modes of the 
individual kinks At the initial impact, some kinetic 
energy is transferred into internal 'shape' modes of the 
kink and antikink. They then separate and propagate 
almost independently, but for initial velocities less than 
Vc they no longer have enough translational energy to 
escape their mutual attraction, and so they return and 
collide a second time. At this point some of the energy 
stored in the shape modes can be returned to the trans- 
lational modes, tipping the energy balance back again 
and allowing the kink and antikink to escape to infin- 
ity, provided that there is an appropriate resonance be- 
tween the interval between the two collisions, and the 
period of the internal modes. More generally, sufficient 
energy might be returned to the translational modes af- 
ter three or more kink-antikink collisions, leading to an 
intricate nested structure of resonance windows, a pic- 



ture which has been confirmed by both numerical and 
analytical studies [9l-[l3j. 

For this mechanism to work, the kink and antikink 
must each support at least one internal vibrational mode, 
within which energy can be stored prior to being trans- 
ferred back to the translational modes after a number 
of bounces. This has led to the commonly- held belief, 
stated for example in 12[ , that the existence of an inter- 
nal kink mode is a necessary condition for the appearance 
of resonance windows. The example of the parametrically 
modified sine-Gordon model lends further support to this 
view: depending on the value of a parameter, kinks and 
antikinks do or do not possess an internal mode; corre- 
spondingly, resonance windows do or do not appear [lo[ . 
Resonances have also been observed in vector soliton col- 
lisio ns fl^ [isl and in the scattering of kinks on impuri- 
ties [16j, but again, the mechanism always relies on the 
presence of a localised internal mode, either of a single 
kink or of an impurity, or both. 

In the present Letter we revisit this question in the 
context of kink-antikink scattering in the cj)^ model . 
The potentials mentioned in [ip| include a general poly- 
nomial form, but the (p^ model has a number of features 
whose relevance to resonant scattering have not been ap- 
preciated to date, even in the recent work 18|. As for the 
(j)'^ model, it is non-integrable, but, as for the modified 
sine-Gordon model in the regime of no resonant scatter- 
ing, a single (j>^ kink does not have an internal oscillatory 
mode 17|. Nevertheless, certain kink-antikink collisions 



exhibit resonant scattering, thus providing a counterex- 
ample to the standard lore. We elucidate the mechanism 
for this process, finding that the underlying reason is a 
reversible transfer of kinetic energy from the kinks to a 
collective bound state trapped by the KK pair. 

The model. The one-dimensional 4>^ theory can be de- 
fined by the rescaled Lagrangian density [ij 
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kinks and antikinks, interpolating between neighbouring 
vacua, can be found from the one-kink solution (px (x) = 
4'{o,i){x) — \/ (l+tanha;)/2 using the discrete symme- 
tries of the model under — J> —0 and/or x — x, so 
= 0(1,0) (a;) = <I>k{~x), 0(o_i)(x) -cpKix) and 
4'(-ifi){x) — —(pKi—x). These all have mass M — 1/4. 

The model also possesses perturbative meson states, 
fluctuations about the vacua (/> = 0, ±1. Fluctuations 
around the static solutions (psix) which interpolate be- 
tween these vacua can be treated by setting 4>(x, t) = 
(t)s{x) + 77(x)e'"*. The linearised field equations are 
—'Hxx + U{x)r] = uj'^r/ where the potential U{x) is 

U{x) = 15(t>t - 120^ + 1- (2) 

For an isolated kink there are no localised solutions to 
this equation beyond the usual translational zero mode, 
reflecting the absence of internal oscillatory modes. The 
states of the continuum spectrum can be written in terms 
of hypergeometric functions ,17] . 

Numerical results. The initial conditions we took cor- 
respond to a widely separated kink-antikink pair propa- 
gating towards a collision point. To find a numerical so- 
lution of the PDF describing the evolution of the system, 
we used a pseudo-spectral method on a grid containing 
2048 nodes with periodic boundary conditions. For the 
time stepping function we used a symplectic (or geomet- 
ric) integrator of 8*'^ order to ensure that the energy is 
conserved. The time and the spatial steps were 6x = 0.25 
and 6t = 0.025. To check numerical stability we repeated 
selected calculations with St = 0.05 and 6t = 0.0125. 

Since the model ([T} contains two distinct classes of 
vacua, and ±1, kink-antikink collisions in the sectors 
built on the (jjv = and the (j)v — vacua have to be 
analysed separately. In the former case the initial config- 
uration, which we denote as KK or (0, 1) -I- (1, 0), can be 
taken as a superposition 0(x) = (fiKix+a) + (pxix—a) — 1 
where a > is the initial (half-)separation parameter, 
which we set equal to 12.5; in the latter case we take 
(l){x) — (j)x(x+a) + (j)K{x—a) and denote this as KK 
or (1,0) -I- (0,1). In neither case are there internal vi- 
brational modes bound to the individual kinks. By the 
standard picture, we would therefore predict that neither 
should exhibit any resonance structure. 

In the case of (0, l)-l-(l, 0) KK collisions, this is indeed 
what we found (Figs.[TJa), (b)). For v <Vcp~ 0.289 the 
pair always become trapped, while for v > Vcp the colli- 
sion yields a reflected pair of solitons together with some 
radiation: (0, 1)-F (1, 0) (0, -1)-F (-1, 0). However the 
(1, 0)-|-(0, 1) KK colhsion reveals a very different picture, 
illustrated in Figs.[T](c), (d) and (e). An intricate pattern 
of escape windows can be seen, up to a critical velocity 
Vcr ~ 0.0457, after which the kinks always have enough 
energy to separate. This is very similar to the behaviour 
of the 0** model, even though there are no internal modes 
of the scattering kinks into which kinetic energy can be 
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Figure 1. {K,K} collisions in the (0,1) + (1,0) ((a), (b)) 
and (1,0) -I- (0,1) sectors ((c), (d), (e)). Plots (a) and (c) 
show the fitted final velocity of the kink as a function of 
initial velocity, with the dotted line indicating the result for 
a purely elastic collision, while (b) and (d) depict the field 
values measured at the collision centre. The final plot, (e), 
shows the times to the first, second and third kink-antikink 
collisions in the (1, 0) + (0, 1) sector. 

transferred. Thus, we have to look for another explana- 
tion of our results. 

The mechanism. The spectrum of linear perturba- 
tions around a single (p^ kink differs from that for the 
0* kink in that the potential U{x) is not symmetrical 
with respect to reflections x ^ —x. As a result the po- 
tential for a kink-antikink pair depends on the order in 
which they appear on the line: a well-separated KK pair 
(Fig-Ilta)) has a raised central plateau, while an equally 
separated KK pair (Fig.^b)) has instead a wide central 
well. If the velocities of the kinks are relatively small, 
we expect - and will justify later - that the adiabatic 
approximation can be used to find the spectrum of small 
fluctuations about these two configurations, using Fq. ([2]) 
for the appropriate choices of (pg . Fig.[3l^b) shows our nu- 
merical results for the (1, 0) -I- (0, 1) KK pair as a func- 
tion of separation parameter a. The two lowest states are 
quasizero modes, rapidly approaching zero from opposite 
sides as the separation grows. On top of these is a tower 
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Figure 2. The potential U(x) (solid lines) for linear per- 
turbations about the composite (0,1) + (1,0) KK (a) and 
(1,0) + (0, 1) KK (b) configurations (dashed lines). 
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Figure 3. (a) Peaks of the Fourier transform of the field at the 
origin in the two-bounce window for v = 0.04548, vs. (b) even 
(solid) and odd (dashed) bound states for the (1,0) + (0, 1) 
configuration, as a function of the half-separation a. 



of meson states with separation-dependent energies. The 
analogous plot for the (0, 1)-|-(1, 0) KK pair, by contrast, 
shows only the two quasizero modes. 

So, a new feature of KK collisions in the 0^ model is 
that although there are no internal modes of the individ- 
ual kinks, energy can be stored after an initial impact in 
the trapped meson states of the composite (1,0) + (0, 1) 
configuration. This opens the possibility that, for colli- 
sions with Vi < Vcr satisfying a suitable resonance condi- 
tion, this energy might be returned to the translational 
modes on a subsequent recoUision, thereby allowing the 
kink and antikink to return to infinity. 

To support this picture, we first assign a 'bounce num- 
ber' to each window, equal to the number of collisions 
that the kink and antikink undergo before their final es- 
cape to infinity. Within each window the bounce number 
is constant. The first two-bounce windows, discernible on 
Fig.[lje) as the intervals of velocity within which the time 
to the third collision diverges, are centred at Vi « 0.0228, 
followed by a 'false' window at Vi w 0.0273, then a third 
at Vi « 0.0303, and so on. The regions near to these 
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Figure 4. (a) Field value at the collision centre as a function 
of time from the first collision for three consecutive two- 
bounce windows, (b) The time between the two collisions in 
the two-bounce windows vs. the number of oscillations. 



windows exhibit nested structures of higher-bounce win- 
dows, similar to those seen in kink-antikink interactions 
model 



the 



11|,|12|. 



Our analysis implies that each two-bounce window 
should be associated with a distinct integer n, the number 
of collective mode oscillations between its two collisions. 
In Fig. lUa) we illustrate this for the three windows just 
discussed. It is easily checked that the number of os- 
cillations of 4>{0, t) between the two collisions increases 
by one from one plot to the next. Furthermore, our nu- 
merical results clearly indicate that these oscillations are 
those of the lowest collective mode of the KK system. 
In the two-bounce windows we considered, we observed 
maximal kink-antikink separations between bounces up 
to R = 2a ~ 12. Fig. ^a) shows the Fourier trans- 
form of the value of the field at the collision centre be- 
tween the times of the first two bounces for Vi = 0.04548. 
The time between these bounces is T « 660, during 
most of which the separation between the kinks is close 
to 12. The Fourier transform shows four peaks at fre- 
quencies 1.045,1.28,1.61,1.92, with the strength of the 
second peak about one order smaller than that of the 
first, and the higher ones strongly suppressed. The posi- 
tions of these peaks relative to the energies of the first 4 
parity-even modes (Fig. |3l[a) ; note that the symmetry of 
the initial condition rules out any coupling to parity-odd 
modes) are consistent with a half-separation between the 
kink and antikink of a « 6.13, in good agreement with 
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the observed value. This also lends support to our use 
of the adiabatic approximation in the analysis of fluctu- 
ations about the kink-antikink configuration. 

In Fig.|3Ub) we show the dependence of the times T{v) 
between the two collisions in the two-bounce windows on 
the number n{v) of oscillations in between these colli- 
sions. The results are a good fit to the linear relation 
uiT = 27rn + 6, where w — 1.0452, the same as the fre- 
quency of the lowest collective mode at 2a « 12 and the 
value found from the Fourier transform at Vi = 0.04548, 
and S = 6.0277. Requiring S to lie between and 27r, 
as here, fixes the number n of collective mode oscilla- 
tions assigned to the three windows shown in Fig. \M.&) 
to be 12, 13 and 14 respectively, while the initial velocity 
Vi — 0.04548 used for the Fourier transform plot is very 
close to the n = 109 two-bounce window. 

A linear fit with this value of w is exactly as predicted 
by a resonance mechanism driven by the first mode, and 
closely resembles similar relations in the 0'* model 0, 
although the detailed mechanism here is different. In 
particular, in resonant cj)'^ scattering, energy is exponen- 
tially localised on the kink and antikink, while in the (p^ 
theory this energy is stored in an extended meson state 
residing in the potential well formed between the kink 
and antikink. This results in a small amount of radiation 
pressure on the kinks after their initial impact, which we 
observed for Vi = 0.05, just larger than Vc, in a small 
(of order 10~^) acceleration of the bouncing kinks away 
from each other, when measured at times t « 300 after 
collision so that the kinks had become well separated. A 
corresponding simulation in 0'* theory showed that any 
such effect is at least 10"* times weaker there. 

Nevertheless, we can pursue the similarity between 
the two models by considering the asymptotic attractive 
forces between kink and antikink in the absence of radi- 
ation pressure, which are F{a) ^ 2e^^ for the KK pair, 
and F{a) ~ 2e~^^ for KK pair, where i? = 2a is the 
kink-antikink separation. The arguments of would 
then predict T{v) oc (w^^ — vf)^°' for impact velocities 
Vi just below Vcr, with a = 0.5. Our fits based on this 
functional form favour a smaller value for a, in the range 
0.39 — 0.45. It would be interesting to see whether a more 
refined treatment, taking radiation pressure into account, 
could account for these discrepancies. 

Perhaps the most intriguing feature of resonant 4>^ 
scattering, though, is the 'missing' window at n = 13. 
For resonant (j>* scattering, two-bounce windows are also 
missing, for n < 3, but once they set in they are found 
for all n, at least up to initial velocities very close to 
Vc 0- By contrast, in (jj^ scattering we found the first 
two-bounce window at n = 12, then a false window at 
n = 13, and then true windows for all higher values of 
n that we examined. Even more remarkably, a similar 
structure is reproduced when looking at the three-bounce 
windows next to a given two-bounce window, and we sus- 
pect that this pattern will continue at all higher levels. It 



is possible that an explanation for this behaviour will be 
found in a careful treatment of the higher modes of the 
bound-state spectrum, but it remains a major challenge 
to convert this idea into a robust set of predictions. 

Conclusions. The intricate structure of (f)^ kink colli- 
sions overturns some previous beliefs about resonant kink 
scattering, and shows that it has wider relevance than 
previously thought. The new mechanism enabling reso- 
nances to occur is the formation of meson bound states in 
the potential well created in the space between the con- 
stituents of a suitably-ordered kink-antikink pair. The 
resulting windows are less regular than for (j>* scattering, 
and exhibit gaps. For large kink-antikink separations we 
also saw an additional long-range interaction through the 
trapped meson field, visible as radiation pressure. 
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